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ABSTRACT 

Planets orbiting post-common envelope binaries provide fundamental information on planet formation and evolution. We searched for 
such planets in NN Ser ab, an eclipsing short-period binary that shows long-term eclipse time variations. Using published, reanalysed, 
and new mid-eclipse times of NN Ser ab obtained between 1988 and 2010, we find excellent agreement with the light-travel-time effect 
produced by two additional bodies superposed on the linear ephemeris of the binary. Our multi-parameter fits accompanied by N-body 
simulations yield a best fit for the objects NN Ser (ab)c and d locked in the 2 : 1 mean motion resonance, with orbital periods _P C a* 
15.5 yrs and Pd -7.7 yrs, masses M c sin ; c as 6.9Mj up and Md sin z'd - 2.2Mj„ p , and eccentricities e c ^0 and -0.20. A secondary x 2 
minimum corresponds to an alternative solution with a period ratio of 5 : 2. We estimate that the progenitor binary consisted of an A star 
with ~ 2 M Q and the present M dwarf secondary at an orbital separation of ~ 1 .5 AU. The survival of two planets through the common- 
envelope phase that created the present white dwarf requires fine tuning between the gravitational force and the drag force experienced 
by them in the expanding envelope. The alternative is a second-generation origin in a circumbinary disk created at the end of this 
phase. In that case, the planets would be extremely young with ages not exceeding the cooling age of the white dwarf of 10 6 yrs. 

Key words. Stars: evolution - Stars: binaries: eclipsing - Stars: individual: NNSer- Stars: cataclysmic variables - Stars: planetary 
systems - Planets and satellites: detection - Planets and satellites: formation 



1. Introduction 

NNSeratQ is a short-period (f 01 b = 3.12hr) eclipsing binary 
at a distance of 500 pc. The detached system contains a hot 
hydrogen-rich white dwarf NNSer a of spectral type DAOl 
and an M4 dwarf star NNSerb with masse s of 0.535 M and 
O.IIIMq, respectively dParsons et al.ll2010al). With an effective 
temperature of 57000 K dHaefner et al. 2004J), the white dwarf 
has a cooling age of only 10 6 yrs (IWoodlll995l) . The present sys- 
tem resulted from a normal binary with a period of ~ 1 year when 
the more massive component evolved to a giant and engulfed the 
orbit of its companion. The subsequent common envelope (CE) 
phase led to the expulsion of the envelope, laying bare the newly 
born white dwarf and substantially shortening the orbital period. 

Some eclipsing post-CE binaries dis play long-term eclips e 
time variations, am ong them V471 Tau (IKa minski et al. 2007), 
QS Vir and NN Ser (Par sons et alj|2010bl and references there- 
in). The latter possesses deep and well-defined eclipses, which 
allow measurement s of the mid-eclipse times to an accuracy o f 
1 00 ms and better dBrinkworth et alj2006t iParsons et all2010bl) . 
The processes advanced to explain them include the long-term 
angular momentum loss by gravitational radiation and magnetic 
braking, possible quasi-periodicities caused, e.g., by Applegate's 

1 On recommendation by the Editor of A&A, we refer to the sys- 
tem as NN Ser, to the binary explicitly as NN Serab, and to the objects 
orbiting the binary as NN Ser (ab)c and NN Ser (ab)d. 



(1992) mechanism, and the strict periodicities produced by apsi- 
dal motion or the presence of a third body in the system. Finding 
the correct interpretation requires measurements of high preci- 
sion and a coordinated effort over a wide range of time scales. 
The existence of a third body orbitin g NN Serab was previ- 
ously considered by lOian et all d2009l) . but the orbital parame- 
ters suggested by th em are incompatible with more recent data 
dParsons et al.l 1201 Obi) . In this Letter, we present an analysis of 
the eclipse time variations of NN Ser ab, based on published 
data, the reanalysis of published data, and new measurements 
obtained over the first half of 2010. 

2. The data 

After their 1 988 discovery of deep eclipses in NNSer, 
lHaefner et al.l d2004l) acquired a series of accurate mid-eclipse 
times in 1989. After a hiatus of ten years, they added a po- 
tentially very accurate trailed CCD imaging observation using 
the ESO VLT. From 2002 on, the Warwick group systemat- 
ically se£ure^_j_Jojal of_22 mid-eclipse times of high preci- 
sion dBrinkworth et al.1 120061: IParsons et all 1201 Obi this work). 
IParsons et al. IHoiOb) list all published mid-eclipse times by 
other authors until the end of 2009. These are included in our 
analysis that weights them by their statistical errors. Since the 
individual Warwick mid-eclipse times between 2002 and 2009 
were separated by about one year, information on eclipse time 



1 



K. Beuermann et al.: Two planets orbiting NN Serpentis 




-2000 



2000 
JD-2450000 



4000 



6000 



1.0 
0.0 
-1.0 



1000 



1.0 
0.0 
-1.0 



2000 



Model 1 



3000 4000 
JD-2450000 

—i i i | i r 



5000 



d 3 



5200 



5250 5300 
JD-2450000 



5350 



Fig. 1. Top: Top: Observed - calculated mid eclipse time dif- 
ferences relative to the best-fit linear ephemeris for Model 1 of a 
single planet orbiting NN Ser ab. Center and bottom: Residuals 
relative to the eccentric-orbit fit for two selected time intervals. 
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Fig. 2. Same as Fig. 1 for Model 2a of two planets orbiting 
NN Ser ab. The contributions of components b and c are indi- 
cated by the dashed and dot-dashed curves, respectively, the 
solid curve shows the combined effect. 



variations on a shorter time scale is lacking. We, therefore, 
organized a collaborative effort of the Gottingen, McDonald, 
and Warwick groups to monitor NN Ser over the first half of 
2010. We used the remotely controlled MONET/North 1.2-m 
telescope at McDonald Observatory via the MONET internet 
remote-observing interface, the McDonald 2.1-m telescope, and 
the ESO 3.5-m NTT. The MONET data were taken in white 
light, the McDonald data with a BG40 filter, and the NTT obser- 
vations were acquired with the ULTRACAM high-speed CCD 
camera equipped with Sloan filters. The mid-eclipse times mea- 
sured in Sloan u', g', and i' are consistent, and we used the g' data 

as the most accurate set for the prese nt purpose. 

The mid-eclipse time derived by Haefner et al. (2004) from 
the trailed VLT image of 1 1 June is the most variant of 

the published eclipse time measurements and was assigned a 
large error of 17 s, although this should be a very precise mea- 
surement, given the very simple form of the eclipses in NN Ser 
and the use of an 8.2m telescope. We reanalysed the image of 
11 June 1999, which started 04:53:05.537 UT with an exposure 
of 1125.7462 s and was taken in good atmospheric conditions. 
The key issue is the conversion of the track from pixel space to 
time. Using two independent method s, we found that the orig- 
inal analysis by Haefne r et alj (120041) was in error and that the 
mid-eclipse time can be determined with an accuracy of 0.20 s 
(cycle E _ 30721). We also r eanalysed the less accurate data of 
iPigulski & Michalskal d2002l) (cycle E = 33233) by including the 
effects of the finite integration times. 

2 http://www.eso.org/public/images/eso9936b/ 



Table 2 lists all previously published, the reanalysed, and the 
new mid-eclipse times shifted to the solar system barycenter and 
corrected for leap seconds. The table also gives the 1 -cr statistical 
errors and the residuals relative to our Model 2a, as shown in 
Fig. 2 and discussed in Sect. 4, below. 

3. The light-travel-time effect in NN Ser 

All measurements of mid-eclipse times of NN Ser ab are dis- 
played in Figs. 1 and 2 as O - C values relative to the model- 
dependent linear ephemerides of the respective fits. Data points 
with errors <1 s and >1 s are shown as green and yellow dots, 
respectively. The eclipse tim e measuremen t s dom inating the fit 
are the 1989 data points of lHaefner et alj d2004 near the ab- 
scissa value JD'= JD-2450000 = -2295, the reevaluated VLT 
point on JD' =1340, the 2002-2009 series of Warwick eclipse 
times since JD' =2411 dParsons et al.ll20i0bh . and the data of 
this work since JD' =5212. In particular, the revised VLT mid- 
eclipse time implies a twofold change in the time derivative of 
O - C and exclu d es the sim ple quadratic ephemeri des used by 
iBrinkworth eTafl d2006h and iParsons et al] d2010bl) . The avail- 
able data do not exclude abrupt period changes or an ultimate 
aperiodicity, but there is no physical process that predicts such 
behavior. We consider a periodic behavior the most promising 
assumption and proceed to explore this possibility. 

Strictly periodic O—C variations may result from apsidal mo- 
tion of the binary orbit or an additional body orbiting the binary. 
Given the parameters of NN Ser ab, classical apsidal motion for 
small eccentricities et,in produces a sinusoidally varying time 
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shift with an amplitude Pb; n ebin/fl' = 3577ebin s (lTodoran|ll972h . 
As a result, <?bi n ~ 0.01 would suffice to produce the observed am- 
plitude. However, the likewise predicted variation of the FWHM 
of the ecli pse and the time shi ft of the secondary eclipse are not 
observed (Parsons et a 1. 2010b, this work). Furthermore, the ob- 
served variation is not sinusoidal and, given an apsidal motion 
constant for the secondary star NNSerb of £22 - 0.11, the pe- 
riod of the apsidal motion would be as short as ~ 0.4 years. Such 
periodicity is not detected (see Fig. 2, bottom panel). 

This leaves us with the third-body hypothesis, at least for the 
major fraction of the observed eclipse time variations. In general, 
it would be possible that different physical processes combine to 
produce the observed signal. We find, however, that a perfect 
fit within the very small statistical errors can be obtained for a 
signal that consists of the periodicities produced by two objects 
orbiting NN Ser ab. Guided by Ockham's razor and the history of 
discoveries in the Solar system, we consider that a fourth body 
in the presence of a third one is a natural assumption. 

4. One-planet and two-planet fits to the data 

Including the light-travel-time effect of the objects NN Ser (ab)c 
and NN Ser (ab)d, the times of mid-eclipse become 

ZA'bink (1 — ^u) 
— ^-sin(u k -zzr k ), (1) 
(l+e k cosu k ) 

k=c,d 

where time is measured from a fiducial mid-eclipse time Tq. A 
linear binary ephemeris is assumed with Pbin the orbital period 
and E the cycle number. The five free parameters for planet k are 
the orbital period i\, the eccentricity e^, the longitude of perias- 
tron tzrbin.k measured from the ascending node in the plane of the 
sky, the time 7\ of periastron passage, and the amplitude of the 
eclipse time variation K\ = flbin.k sin i^/c, with flbin.k the semi- 
major axis of the orbit of the center of mass of the binary about 
the common center of mass of the system, i\ the inclination, and 
c the speed of light. In the denominator, is the true anomaly, 
which progresses through 2n over the orbital period i\- 

We explored the multi-dimensional x 2 space of the two- 
planet model, using the Levenberg-Marquardt routine imple- 
mented in IDL and an independent code. The search showed that 
compensation effects render some parameters ill defined. This 
uncertainty results, in particular, from the long hiatus bet ween 
the accurate measurements of 1989 (Hae fner et alJ [2004) and 
1999 (VLT, this work). We selected the best model, therefore, 
by imposing the additional requirement that the derived orbits 
be secularly stable. We investigated all solutions permitted by 
the data with numerical N-body simulations with a variable time 
step Runge-Kutta integrator, following the orbits over 10 5 yrs, 
and find that only a narrow range in parameter space corresponds 
to stable solutions. In what follows, we consider the one-planet 
and the two-planet models in turn. 

Model 1 with seven free parameters describes a single planet 
with eccentricity e. The fit requires e ~Z 0.60 and is bad for any 
value of e, with a reduced^ >23.3 (x 2 = 1052 for 45 degrees 
of freedom). The top panel of Fig 1 shows the case <? = 0.65. The 
residuals based on the statistical errors of the data points (cen- 
ter panel) reach 23 standard deviations and indicate that there 
is an additional modulation at about half the orbital period. The 
residuals of the 2010 data (bottom panel) demonstrate the lack 
of O - C fluctuations on a short time scale. 

Model 2 for two planets requires some restriction in param- 
eters, because the grid search yields good fits for a range of ec- 
centricities of the outer planet <? c , including zero, and for a period 



ratio r p = PJP d = 1.90 ± 0.30 or r p = 2.50 + 0.15 (1-cr errors), 
with the former slightly preferred. The dichotomy in r p arises 
from the uncertain phasing of the singular 1989 point relative to 
the train of the 1999-2010 data. Further minima at still larger 
r p do not exist. Only a small fraction of the parameter space al- 
lowed by the fits corresponds to secularly stable orbits, however. 
Near r p ^ 2, orbits with <? c > 0. 1 tend to be unstable, while the sta- 
bility region is broad in the remaining parameters for <? c w 0.02. 
Furthermore, all solutions with r p $ 1.9 are unstable, with only 
some solutions stable at r p = 1.9. The solutions near r p = 2.5 are 
more generally stable. We consider Models 2a and 2b, represent- 
ing the cases of r p 2.0 and 2.5, respectively, both with e c = 0. 
Model 2a provides the slightly better fit and is shown in Fig. 2. 
It yields K c = 27.4 s, K d = 5.7 s, P c = 15.5 yrs, P A = 7.75 yrs, and 
e A = 0.20 with xl = 0.78 (x 2 = 32.9 for 42 d.o.f .)• Periastron pas- 
sage of NNSer(ab)d occurred last on JD'^ 45 15. At that time 
NNSer(ab)c was at longitude 213°. For the low value of e c = 
0.03, a shallow minimum of x 2 is attained for aligned apses. 
From the present data, we cannot infer the true value of r p with 
certainty, but it is intriguing that objects c and d may be locked in 
either the 2 : 1 resonance, found also in other planetary systems, 
or the 5 : 2 resonance. The parameters for Models 2a and 2b are 
listed in Table 1, together with their 1-cr errors. A simpler model 
with two circular orbits reaches only xl = 1 -96 (x 2 = 86.2 for 44 
d.o.f.) at r p =2.46 and can be excluded. 

Using Model 2a as input to our N-body simulations, we find 
that <? c and e& oscillate around 0.02 and 0.22 with amplitudes 
of 0.02 and 0.05, respectively. The difference Act of the perias- 
tron longitudes circulates on a time scale of 400 yrs. The periods 
perform small-amplitude anti-phased oscillations, which cause 
r p to oscillate between 1.9 and 2.2. Even if the two planets are 
secularly locked in the 2 : 1 mean motion resonance, therefore, 
the observed period ratio at any given time may deviate slightly 
from its nominal value. 

For Model 2a, the best-fit binary ephemeris is T — BJED 
2,447344.524425(40) + 0.1300801419(10)£, where the errors 
refer to the last digits. Adding a quadratic term BE 1 to the 
ephemeris does not improve the two-planet fit and yields a 1 -<x 
limit of I B I < 1 .5 10 -13 days, leaving room for a period change 
by gravitational radiation or a long-term ac tivity-related effect 
(Bri nkworth et al.ll2006t|Parsons et alj|2010bl) . 

5. Discussion 

The large amplitude of the O - C eclipse time variations in 
NN Ser can only be explained by a third body in the system, 
while the still substantial residuals from a single-planet fit could, 
in principle, have a different origin from that of a fourth body. 
The two-planet model, however, possesses the beauty of simplic- 
ity, and the fact that the residuals for the entire data set vanish si- 
multaneously imposes tight restrictions on any other mechanism. 
In particular, the lack of short-term variability of the residuals in 
the first half of 2010 argues against any process that acts on a 
short time scale or leads to erratic eclipse time variations. Hence, 
there is strong evidence for two planets orbiting NN Ser ab. 

With masses M c sin/ C 6Mj up and Mdsin/d - 2Mj up , 
NN Ser (ab)c and NN Ser (ab)d both qualify as giant planets for 
all inclinations z' c >28° and i& >9°, respectively. The probable 
detection of resonant motion with a period ratio of either 2 : 1 or 
5 : 2 is a major bonus, which adds to the credence of the two- 
planet model. It is the second planetar y system found by eclipse 
timing, after HW Vir dLee et alJl2009h . 

Given a pair of planets orbiting a post-CE binary, two for- 
mation scenarios are possible. They could either be old first- 



3 



K. Beuermann et al.: Two planets orbiting NN Serpentis 
Table 1. Parameters of the models fitted to the measured mid-eclipse time variations of NN Ser, where = indicates a fixed parameter. 



Model 


Planets 


Number 


Pc 


Pd 


PJPi 






a c 


ad 






M c sin i c 


M d sin ; d 


x 2 


xl 






free par. 


(yrs) 


(yrs) 








(AU) 


(AU) 


O 


O 


(M, up ) 


(M, up ) 






1 


1 


2 + 5 


22.60 






>0.65 




6.91 




8.0 




8.36 




1052.3 


23.38 


2a 


2 


2 + 8 


15.50 


7.75 


2.00 


= 0.0 


0.20 


5.38 


3.39 




74 


6.91 


2.28 


32.9 


0.78 








±0.45 


±0.35 


±0.15 




±0.02 


±0.20 


±0.10 




±4 


±0.54 


±0.38 






2b 


2 


2 + 8 


16.73 


6.69 


2.50 


= 0.0 


0.23 


5.66 


3.07 




73 


5.92 


1.60 


33.8 


0.80 








±0.26 


±0.40 


±0.15 




±0.04 


±0.06 


±0.13 




±7 


±0.40 


±0.27 







generation planets that formed in a circumbinary protoplanetary 
disk or they could be young second-generation p lanets form ed % 
10 6 yrs ago in a disk that resulted from the CE (lPeretsll20 1 0b . To 
evaluate both scenarios, we have reconstructed t he CE evolution 
of NN Ser ab using the improved algorithm by IZorotovic et al.l 
d2010l) . who constrain the CE efficiency to a range a = 0.2 - 0.3. 
Possible solutions for the progenitor binary of NN Ser ab are 
not very sensitive to a: for a = 0.25, the progenitor was a gi- 
ant of 2.08 M Q and radius 194 R Q with the present secondary 
star at a separation of 1 .44 AU. When the CE engulfed the sec- 
ondary star, dynamic friction caused the latter to spiral in rapidly, 
thereby dramatically decreasing the binary separation to the cur- 
rent 0.0043 AU. Stability arguments imply that any planet from 
the pre-CE phas e must have formed with semi-major axes ex- 
ceeding 3.5 AU (Holman & Wiegerj [T999t) . With three quarters 
of the central mass expelled in the CE event, pre-existing plan- 
ets would move outward or may even be lost from the system. 
However, given a sufficiently dense and slowly expanding CE, 
the dynamical force experienced by them m ay have ultimately 
moved them inward dAlexander et alj[l976l) . Since the drag pri- 
marily affects the more massive and more slowly moving outer 
planet, such a scenario could lead to resonant orbits, so a first- 
generation origin appears possible. 

The alternative post-CE origin in a second-generation of 
planet formation is also possible, since the formation of cir- 
cumbinary disks is a common phenomenon among post-AGB 
binary stars and the concentration of a slow, dusty wind to the 
orbital plan e of the binary is thought to favor th e formation of 
planets (e.g. lvan Winckel et al.ll2009t lPeretsll20 1 0l) . In particular 
the tiny separation of the present binary poses no problem for 
stable orbits of second-generation planets even at significantly 
shorter distances than the inner planet that we have detected 
dHolman & Wieger3ll999l) . A particularly intriguing aspect of 
a second-generation origin of the planets in NN Ser would be 
their extreme youth, eq ual to or less than the 10 6 yrs cooling 
age of the white dwarf (IWoodlll995h . This feature would dis- 
tinguish them from all known exoplanets and may ultimately 
lead to their direct detection. While we cannot presently prove 
a second-generation origin for these planets, modeling the CE 
event may allow us to distinguish between the two scenarios. 
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Table 2. Previously published, reanalysed, and new mid-eclipse times of the white dwarf in NN Ser with residuals for the light- 
travel-time effect produced by the two planets of Model 2a. The published mid-eclipse times have been converted to BJD(TT) if not 
yet on this time standard. 



E 


BJD(TT) 


Error 


Residual 


Error 


Residual 


References 


Comment 




JD2400000+ 


(days) 


(days) 


(s) 


(s) 









47344.5246635 


0.0003500 


0.0000290 


30.00 


2.51 


(1) 


Reanalysed 


2760 


47703.5457436 


0.0000020 


0.0000012 


0.17 


0.10 


(2) 




2761 


47703.6758326 


0.0000060 


0.0000101 


0.52 


0.87 


(2) 




2769 


47704.7164596 


0.0000030 


-0.0000038 


0.26 


-0.33 


(3) 




2776 


47705.6270226 


0.0000030 


-0.0000016 


0.26 


-0.14 


(3) 




2777 


47705.7571046 


0.0000070 


0.0000003 


0.60 


0.03 


(3) 


Corrected 


2831 


47712.7815836 


0.0001500 


0.0001534 


12.96 


13.25 


(2) 




2839 


47713.8222336 


0.0001500 


0.0001625 


12.96 


14.04 


(2) 




7360 


48301.9141954 


0.0001500 


-0.0000627 


12.96 


-5.42 


(2) 




28152 


51006.5405495 


0.0002000 


0.0000605 


17.28 


5.23 


(2) 




30721 


51340.7165402 


0.0000023 


-0.0000004 


0.20 


-0.03 


(2) 


Reanalysed 


33233 


51667.4780058 


0.0000960 


0.0000041 


8.29 


0.35 


(4) 


Reanalysed 


38960 


52412.4470566 


0.0000006 


-0.0000006 


0.05 


-0.05 


(5) 




38961 


52412.5771382 


0.0000005 


0.0000008 


0.04 


0.07 


(5) 




38968 


52413.4876977 


0.0000009 


-0.0000006 


0.08 


-0.05 


(5) 


Corrected 


38976 


52414.5283389 


0.0000007 


-0.0000004 


0.06 


-0.03 


(5) 




38984 


52415.5689804 


0.0000007 


0.0000000 


0.06 


0.00 


(5) 




41782 


52779.5331703 


0.0000015 


0.0000001 


0.13 


0.01 


(5) 




41798 


52781.6144523 


0.0000007 


0.0000002 


0.06 


0.02 


(5) 




41806 


52782.6550927 


0.0000008 


-0.0000004 


0.07 


-0.03 


(5) 




41820 


52784.4762150 


0.0000008 


0.0000003 


0.07 


0.03 


(5) 




44472 


53129.4486808 


0.0000040 


0.0000008 


0.35 


0.07 


(5) 




44473 


53129.5787632 


0.0000028 


0.0000031 


0.24 


0.27 


(5) 




44474 


53129.7088370 


0.0000017 


-0.0000032 


0.15 


-0.28 


(5) 




44480 


53130.4893234 


0.0000030 


0.0000025 


0.26 


0.22 


(5) 




49662 


53804.5644567 


0.0000025 


0.0000001 


0.22 


0.01 


(5) 




49663 


53804.6945350 


0.0000012 


-0.0000017 


0.10 


-0.15 


(5) 




49671 


53805.7351781 


0.0000006 


0.0000005 


0.05 


0.04 


(5) 




53230 


54268.6903114 


0.0000006 


0.0000008 


0.05 


0.07 


(5) 




53237 


54269.6008713 


0.0000002 


-0.0000001 


0.02 


-0.01 


(5) 




56442 


54686.5076279 


0.0000009 


-0.0000001 


0.08 


-0.01 


(5) 




58638 


54972.1634971 


0.0000800 


-0.0000380 


6.91 


-3.28 


(6) 




58645 


54973.0740553 


0.0001000 


-0.0000406 


8.64 


-3.51 


(6) 




58684 


54978.1471791 


0.0001200 


-0.0000408 


10.37 


-3.53 


(6) 




58745 


54986.0820789 


0.0001200 


-0.0000274 


10.37 


-2.37 


(6) 




58753 


54987.1228359 


0.0001300 


0.0000887 


11.23 


7.66 


(6) 




58796 


54992.7161925 


0.0000015 


0.0000008 


0.13 


0.07 


(6) 




60489 


55212.9418187 


0.0000069 


0.0000027 


0.60 


0.23 


(7,8) 




60505 


55215.0230961 


0.0000066 


-0.0000017 


0.57 


-0.15 


(7,8) 




60528 


55218.0149380 


0.0000043 


-0.0000024 


0.37 


-0.21 


(7,8) 




60735 


55244.9415254 


0.0000029 


0.0000012 


0.25 


0.10 


(7,8) 




60743 


55245.9821654 


0.0000032 


0.0000003 


0.28 


0.03 


(7,8) 




60751 


55247 0228063 


0000034 


0000002 


0.29 


0.02 






60774 


55250.0146469 


0.0000034 


-0.0000018 


0.29 


-0.16 


(7,8) 




60927 


55269.9169047 


0.0000014 


-0.0000018 


0.12 


-0.16 


(7,9) 




60950 


55272.9087487 


0.0000013 


-0.0000005 


0.11 


-0.04 


(7,9) 




61219 


55307.9003015 


0.0000010 


0.0000005 


0.09 


0.04 


(7,10) 




61426 


55334.8268834 


0.0000018 


-0.0000025 


0.16 


-0.22 


(7,9) 




61440 


55336.6480059 


0.0000018 


-0.0000017 


0.16 


-0.15 


(7,9) 




61441 


55336.7780894 


0.0000015 


0.0000017 


0.13 


0.15 


(7,9) 




61564 


55352.7779443 


0.0000016 


0.0000017 


0.14 


0.15 


(7,9) 




61579 


55354.7291448 


0.0000009 


0.0000004 


0.08 


0.03 


(7,10) 





(1) Haefner, R., 1989, ESO Msngr, 55, 61, reanalysed using up-to-date eclipse profile; (2) Haefner et al. (2004), misprint for E=2777 corrected, 
VLT trailed imaging observation (E=30721) reanalysed using the original data; (3) Wood, J. H. & Marsh, T. R., 1991, ApJ, 381, 551; (4) Pigulski 
& Michalska (2002), reanalysed using the original data; (5) Parsons et al. (2010b), timing for E=38968 corrected for misprint; (6) Qian et al. 
(2009); (7) This work, (8) MONET/North 1.2-m white light photometry, (9) McDonald 2.1-m photometry with Schott BG40 filter, (10) ESO 
NTT 3.5-m ULTRACAM Sloan g' photometry. 
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